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ABSTRACT 

We have used a coupled dynamical and chemical model to examine the chemical 
changes induced by the passage of an interstellar shock in well shielded regions. Using 
this model we demonstrate that the formation of H2O in a shock will be followed in 
the post-shock phase by depletion of the water molecules onto the grain surfaces. To 
attempt to discriminate between the creation of ices behind shocks and their production 
by means of grain surface chemistry, we examine the deuterium chemistry of water 
before, during, and after a shock. We show that chemical evolution in the post-shock 
gas can account for both the deuterium fractionation and the abundance of CO2 
relative to H2O observed in interstellar and cometary ices. Given the pervasiveness of 
shocks and turbulent motions within molecular clouds, the model presented here offers 
an alternate theory to grain surface chemistry for the creation of ices in the interstellar 
medium, ices that may ultimately be incorporated into comets. 

Subject headings: comets: general — ISM: abundances — ISMxlouds — ISM:molecules 
- molecular processes — shock waves 
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1. Introduction 

There is mounting evidence that some interstellar material survives inside objects in the 
solar system. Such evidence is found in the compositional similarity between cometary ices and 
inter-planetary grains with those found in the interstellar medium (Mumma 1997). That cometary 
ices and organics may have supplied the molecular seeds for prebiotic life on the earth, or even 
the oceans, leads to the intriguing possibility that the molecules essential for life might have 
their origin in interstellar space (e.g. Chyba & Sagan 1997). As such it is important to gain an 
understanding of the processes that lead to the formation of ices in the interstellar medium (ISM). 

By far, the most abundant constituent of ices in the interstellar medium, and also in comets, 
is H2O. Because the low temperature ion-molecule chemistry that is active in dense regions of the 
ISM is incapable of reproducing the observed H2O ice observations (Jones Sz Williams 1984), it 
has been proposed that water-ice is formed through the hydrogenation of oxygen atoms on the 
surfaces of cold (Td ~ 10 K) dust grains (c.f. Tielens & Hagen 1982). Until recently there have 
been few alternatives to grain surface chemistry to account for the abundance of ices. However, 
we have found that water-ice mantles form quite naturally in a layer of gas rich in water vapor 
that has been processed by a shock with a velocity in excess of 10 kins -1 (Bergin, Melnick, & 
Neufeld 1998; hereafter BMN). The prediction that large quantities of water vapor are produced 
in shocked gas has recently gained additional support through the direct observations of strong 
emission from gaseous water toward Orion BN-KL (Harwit et al. 1998) using the Infrared Space 
Observatory. 

There are several observational constraints that might discriminate between water-ice mantles 
created behind shocks and those formed by chemical processes in grain mantles. In this Letter we 
examine whether post-shock chemistry can account for the observed abundances of deuterated 
water (HDO) and CO2 in the ISM and in comets. 

2. Model 

We use the 3-stage model described in BMN to examine the evolution of chemical abundances 
in pre-shocked (stage 1), shocked (stage 2), and post-shock gas (stage 3). The chemical evolution 
in each stage is treated independently, except that the final composition of the preceding stage 
is used as the initial composition of the next. The first stage is assumed to be initially atomic 
in composition and has physical conditions appropriate for quiescent gas in molecular cores 
(T = 10 — 30 K; n# 2 = 10 4 — 10 6 cm -3 ). The stage 1 pre-shock chemistry evolves until t = 10 6 yr 
when the gas is assumed to be shocked. The stage 2 shock dynamics are modeled as an increase in 
the gas temperature which scales with the shock velocity — peak gas temperatures range between 
400 and ~2000 K for shock velocities, v s , between 10 and 30 kms -1 , respectively (Kaufman & 
Neufeld 1996). The high-temperature chemistry in the shock is allowed to evolve for a cooling 
timescale (~ 100 yr for a 30 kms -1 shock; see BMN), whereupon the third stage commences with 



-3- 



a return to quiescent temperature conditions. 

We use the UMIST RATE95 database (Millar, Farquhar, & Willacy 1997) for all computations 
of non-deuterated molecules. For the deuterium chemistry we have created a smaller network 
containing all the important species and reactions that lead to the formation and destruction of 
HDO and OD. The database was created using reactions taken from the literature (Croswell & 
Dalgarno 1985; Millar, Bennet, & Herbst 1989; Pineau des Forets, Roueff, & Flower 1989; Rodgers 
& Millar 1996). The principle reaction in the deuterium chemistry is + HD «-> H2D + + H2, 
for which we have used the rate coefficients in Caselli et al. (1998). We have adopted the H3" 
electron recombination rate and branching ratios given in Sundstrom et al. (1994) and Datz et al. 

(1995) , respectively. The H2D + recombination rate and branching ratios are from Larsson et al. 

(1996) . The deuterium chemical network was tested against the larger network of Millar, Bennet, 
& Herbst (1989) and the degree of deuterium fractionation in H2D + , OD, and HDO was found to 
be in excellent agreement for the temperature ranges considered here. As a corollary to the high 
temperature reactions which form H2O in the shocked gas (see Kaufman & Neufeld 1996), we 
have included similar reactions which form HDO. These reactions and rate coefficients are listed 
in Table 1. 

We have used the gas-grain adaptation of the UMIST network discussed in Bergin k, Langer 

(1997) . For the first and second stages we assume that molecules are depleting onto bare silicate 
grains, while in the third stage the molecules deplete onto a grain mantle dominated by solid 
H2O. This will increase the binding potential for all species by a factor of 1.47. We have used 
the measured binding energies of CO, CO2, and H2O to H2O given in Sandford & Allamandola 
(1990). The dust temperature, which is critically important for the rate of thermal evaporation, 
is assumed to be equivalent to the gas temperature, except in the shock (stage 2) where we 
artificially raise the temperature to 200 K to account for the removal of the ice mantle due to 
sputtering or grain-grain collisions. 

3. Results 

Figure 1 presents the time evolution of abundances for the 3-stage gas-grain chemical model. 
The peak gas-phase H2O abundance is low ([H20]/[H2] = 3 x 10~ 7 ) in the pre-shock gas, rises to 
~ 10~ 4 following the passage of a shock, and depletes onto the grain surface at t ~ 10 5 yr (for 
nn 2 = 10 5 cm~ 3 ) in the post-shock gas. In the shock stage we assume v s = 20 kms -1 , which 
will heat the gas to ~1000 K; similar results would be found for any shock velocity between ~10 
and 40 kms -1 . In the pre-shock quiescent stage we find significant levels of D-fractionation: at 
t = 10 6 yr we find [HDO]/[H 2 0] ~ 10~ 3 in both the gas and solid phases. Thus [HDO]/[H 2 0] 
> [HD]/[H2] = 2.8 x 10 -5 , which is the result of low temperatures and low electron abundances 
favoring production of H2D + relative to H3". These enhancements are mirrored in the daughter 
products of H2D + , such as HDO. In the shock itself HDO is released from the grain surface leading 
to an increase in its gas phase abundance. Thus, in stage 2 the HDO abundance shows little 
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change because the fractionation via the ion-molecule chemistry is halted at higher temperatures. 
There is some production via the high-T reactions listed in Table 1; however, these reactions do 
not increase the fractionation significantly. Therefore, the rapid hydrogenation of oxygen is not 
followed by a similar increase in the HDO abundance and in the shock [HDO]/[H20] = 5 x 1CP 5 . 
Pineau des Forets et al. (1989) examined the D-chemistry in shocks and found little change in 
the water D/H ratio. As noted by the authors, their work did not examine the high temperatures 
required to produce abundant H2O and therefore did not probe this change. 

In the post-shock gas (stage 3 in Figure 1) the abundance of gaseous HDO is approximately 
constant with a small, but increasing, abundance frozen on the grain surfaces. At later times, 
t ~ 4 x 10 4 yr, H2O begins to deplete onto the grain surfaces, removing the primary destruction 
pathway for and H2D + As the abundance of H^D" 1 " ions increase, they react directly with 
water forming HDO (along with OD and O), and further increase the fractionation. 

As shown the top panels of Figure 1, the abundance of DCO + is quite high with 
[DCO+]/[HCO+] = 0.002, but in the shock the abundances of both DCO+ and HCO+ decline 
(due to reactions with H2O) and fractionation is reversed such that [DCO + ]/[HCO + ] ~ 2 x 10~ 5 
(see also Pineau des Forets et al. 1989). Thus, DCO + should be an excellent tracer of quiescent 
gas in star forming regions, as chemical theory predicts it will be destroyed in an outflow. 

In the first two stages active CO sublimation at = 30 K suppresses the net CO depletion 
onto bare silicate grains. In stage 3, at ~ 10 5 yr, the gaseous CO abundance declines because we 
have assumed that in the third stage the grains are a priori coated by water-ice which results in a 
greater binding energy, resulting in an exponentially decreasing sublimation rate compared to bare 
silicates (see §2). This leads a solid CO abundance given by [CO] ffr /[H2] = 8 x 10~ 5 at t = 10 7 yr. 
Second, on similar timescales the abundance of OH is enhanced due to the destruction of water 
molecules. The reaction of OH with CO will form CO2, which readily depletes onto grain surfaces. 



4. Discussion 

4.1. Comparison with Observations 

Observations of ices in the ISM and in comets have yielded limits on the level of deuterium 
fractionation in water along with the quantity of CO2 frozen in the water matrix (see Whittet et 
al. 1998). One question is whether our model can simultaneously reproduce both the observed 
[C02]/[H20] and [HDO]/[H20] ratios. Equilibrium models of diffusion-limited grain surface 
chemistry (i.e. chemistry that is limited by the depletion rate onto the surface) can reproduce 
high D-fractionation of water (Tielens 1983) and produce [C0 2 ]/[H 2 0] ~ 0.01 - 10% (Tielens & 
Hagen 1982; Shalabiea, Caselli, & Herbst 1998). 

In Figure 2 we present the comparison of observations with our model predictions as a function 
of time. To illustrate the sensitivity of these results to our assumed gas and dust temperature of 
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30 K, results for higher (40 K) and lower (25 K) gas and dust temperatures are also shown. The 
small circles (with error bars) depict HDO and CO2 measurements for Comets Halley, Hyakutake, 
and Hale-Bopp, the only comets for which all of these quantities have been determined. [] In 
the ISM there are no similar complementary measurements; HDO and H2O have been observed 
in the gas phase in molecular hot cores, while CO2 ices have only been observed along single 
lines of sight towards stars behind molecular clouds or embedded within them. We have therefore 
represented these two independent measurements for the ISM as a dotted box in Figure 2. In hot 
cores, the gas-phase D-fractionation is believed to trace the solid (D/H) ratio because the levels 
of deuterium fractionation are larger than expected for pure gas-phase chemistry evolving at the 
observed temperatures (T ~ 100 — 200 K). Thus, the observed HDO and H2O are proposed to be 
"fossilized" remnants from a previous low temperature phase, remnants that were frozen on the 
grains and evaporate in the hot core (eg. Gensheimer et al. 1996). 

The contours of time-dependent abundance ratios for the range of dust temperatures provided 
in Figure 2 nicely encompass the observations. For this model to be applicable to interstellar ices, 
the initial gas clouds that collapsed and ultimately formed the molecular hot cores must have 
evolved at T ~ 25 — 40 K to account for the observed fractionation. For both CO2 and HDO 
production and where CO2 ices have been observed, our model also suggests that the quiescent 
clouds that formed hot cores would have undergone a shock(s) within the past 10 6 — 10 yr. This 
timescale is consistent with the range of inferred shock timescales in molecular clouds (Reipurth 
et al. 1998; BMN). 

Over the past twenty years or more, there has been considerable debate about whether comets 
are composed of relative pristine interstellar material or whether significant processing has taken 
place within the proto-solar nebula (or indeed whether comets represent a mixture of pristine 
and processed material). The striking similarity between the composition of comets and that of 
interstellar ices, and in particular the very similar level of deuterium fractionation, might suggest 
that the degree of processing in the proto-solar nebula is relatively small. That the composition of 
interstellar ices was preserved during their passage into the outer parts of the proto-solar nebula 
is perhaps unsurprising given the relatively benign physical conditions; in particular, the typical 
shock velocity during accretion onto the trans-Neptunian region (where it is believed that the 
comets originally formed; e.g. Safronov 1969) was only ~ 3kms _1 (e.g. Neufeld & Hollenbach 
1994), too small to result in significant water production in the gas phase. 



Note that while the three comets have similar [HDOJ/pHbO] ratios, Comet Hale-Bopp shows a larger [CC^l/pHhO] 
ratio than either Halley or Hyakutake. These differences are significant relative to the observational uncertainties, 
but may be an artifact of the various Sun-comet distances at which the measurements were made: ISO observations of 
CO2 in Hale-Bopp were obtained when the comet was 2.9 AU from the Sun while the other cometary measurements 
were made at a Sun-comet distance of ~ 1 AU. The lower binding energy of CO2, compared to H2O, could therefore 
lead to differences in the out-gassing of these species at larger Sun-comet distances. Interestingly, ISO observations 
of the short period Comet 103P/Hartley 2 at 1 AU from the Sun also indicate high CO2 abundances, but this comet 
does not have a measured D/H ratio (Crovisier 1998). 
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4.2. Sensitivity to Physical Conditions 

Although our model assumes n# 2 = 10 5 cm -3 , the results are similar at both higher and 
lower densities because most chemical processes show the same dependence on density. We adopt 
Ay = 10 mag in our models, which are therefore only applicable in regions where the ultraviolet 
field is heavily attenuated. Variations in the dust temperature can affect our results, and our 
model can only create ice mantles at dust temperatures greater than 25 K. Below 25 K, O2 
remains frozen on the grain surface in the pre-shock quiescent stage leaving little free gas-phase 
O for shock chemistry to form H2O in abundance. Under these conditions a problem arises in 
accounting for the high abundance of water-ice (~ 10~ 4 ; Schutte 1998) observed in the cold 
foreground gas towards Elias 16. However, this limitation can be mitigated by adopting a lower 
density n# 2 < 10 4 cm~ 3 , as the reduced depletion rate results in a smaller and less constraining 
O2 abundance on grains. 

One unresolved issue concerns the large abundance of CO ice predicted in our model, which 
is quite high (~ 70% relative to H2O) and is larger than observed in the ISM or in comets where 
abundances are generally < 20% (e.g. Mumma 1997). However, if we consider a single CO 
molecule approaching a surface with nearly equal amounts of CO and H2O, it is difficult to imagine 
that CO will always be sharing a (stronger) physisorbed bond with H2O rather than a (weaker) 
bond to a CO molecule. In laboratory experiments where CO and H2O (with CO/H20= 0.5) are 
co-deposited onto metal substrate at 10 K, the evaporation rate of amorphous CO is dominated by 
the vapor pressure of pure CO and is not affected by the presence of H2O (Kouchi 1990). Although 
the sublimation of pure CO is not included in our model (we include sublimation of CO embedded 
in a water mantle), in reality some co-deposited CO will evaporate off the grain mantle, provided 
that the dust temperature exceeds 17 K, the sublimation temperature of pure CO ice. The large 
abundance of frozen CO predicted by our model (under the assumption that CO sublimation is 
negligible) may therefore be unrealistic and should be regarded as an upper limit. 

5. Conclusion 

In this Letter we have demonstrated that in well shielded regions the gas phase chemistry 
after the passage of a shock is capable of reproducing the observed abundance of frozen water, 
the HDO/H2O ratio, and the abundance of CO2 ice. Thus, the creation of grain mantles behind 
shock waves is a viable alternate theory to the production of ices through grain surface chemistry. 
We note that this mechanism does not preclude grain surface chemistry and indeed it may act 
as a supplement. A future extension of our model will be to determine whether shock models 
can account for the enhanced abundances of methanol observed in shocked regions relative to 
surrounding gas (e.g. Bachiller et al. 1995). Such an extension of our model might explain the 
formation of methanol ices - an observed constituent of interstellar grain mantles - as arising from 
the "freeze-out" of CH3OH gas produced by gas phase chemistry in the warm shocked gas. Other 
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tests could involve an examination of other proposed outflow chemistry tracers (e.g. SO) and an 
expanded examination of the deuterium chemistry, including species such as DCN, HDCO, and 
D 2 CO. 

It has previously been suggested that the high degree of deuterium fractionation observed in 
cold cloud cores can explain the observed deuterium fractionation in cometary ices (e.g. Geiss 
&: Reeves 1981). However, our results indicate that shock waves are a critical ingredient in 
the production of ices from gas phase chemistry. Ion-molecule chemistry in quiescent gas can 
reproduce the degree of deuterium fractionation observed in molecular hot cores; however, it 
cannot reproduce the observed abundance of water-ice (Jones & Williams 1984; Bergin, Langer, 
& Goldsmith 1995). In our model, the abundance of water-ice in the pre-shock gas at 10 6 yr 
is < 10 -5 , much less than the value estimated along the line of sight towards Elias 16 (~ 10~ 4 ; 
Schutte 1998). At n# 2 = 10 5 cm~ 3 , the gas-phase chemistry could, with time, eventually approach 
the measured abundance. However, the molecular ice absorption observations are taken along a 
lines of sight that likely have lower densities, which makes reproducing the observed abundances 
more difficult. Thus, either grain surface chemistry, creation of ice mantles behind shock waves, or 
other unknown theories, are necessary to account for the presence of molecular ices. 
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temperature reactions for deuterium-bearing molecules. E.A.B. is also grateful to J. Crovisier 
for discussions of cometary ice abundances. E.A.B. acknowledges support from NASA's SWAS 
contract to the Smithsonian Institution (NAS5-30702); and D.A.N, acknowledges the support of 
the Smithsonian subcontract SV-62005 from the SWAS program. 

REFERENCES 

Bachiller, R. & Perez Gutierrez, M. 1997, ApJ, 487L, 93 
Bachiller, R., Liechti, S., Walmsley, C. M., & Colomer, F. 1995, A&A, L51 
Balsiger, H., Altwegg, K., & Geiss, J. 1995, J. Geophys. Res., 100, 5827 
Bergin, E.A., Melnick, G.J., & Neufeld, D.A. 1998, ApJ, 499, 777 
Bergin, E. A., & Langer, W. D. 1997, 486, 316 

Bergin, E. A., Langer, W. D., & Goldsmith, P. F. 1995, ApJ, 441, 222 
Bockelee-Morvan, D., et al. Icarus, 133, 147 

Caselli, P., Walmsley, CM., Terzieva, R., & Herbst, E. 1998, ApJ, 499, 234 

Chyba, CF. & Sagan, C. 1997, in "Comets and the Origin and Evolution of Life", eds. P.J. 
Thomas, C.J. Chyba, & CP McKay (New York; Springer), 147 



-8- 



Croswell, K. & Dalgarno, A. 1985, ApJ, 289, 618 
Crovisier, J. 1998, priv. comm. 
Crovisier, J. et al. 1997, Science, 275, 1904 
Datz, S., et al. 1995, Phys. Rev. Lett., 74, 896 

Eberhardt, P., Reber, D., Kranowsky, D., & Hodges, R.R. 1995, A&A, 302, 301 
Geiss, J. & Reeves, H. 1981, A&A, 93, 189 

Gensheimer, P.D., Mauersberger, R., & Wilson, T.L. 1996, A&A, 314, 281 
Harwit, M., Neufeld, D.A., Melnick, G.J., & Kaufman, M.J. 1998, ApJ, 497, L105 
Joseph, T., Truhlar, D.G., & Barrett, B.C. 1988, J. Chem. Phys., 88, 6982 
Jones, A.P & Williams, D.A. 1984, MNRAS, 209, 955 
Kaufmann, M. J., & Neufeld, D. A. 1996, 456, 611 

Kranowsky, D. 1991, in Comets in the Post-Halley Era, eds. R.L. Newburn, Jr., M. Neugebauer, 
& J. Rahe (Dordrecht: Kluwer), 831 

Kouchi, A. 1990, J. Crystal Growth, 99, 1220 

Larsson, M., et al. 1996, A&A, 309, LI 

McPhate, J.B., Feldman, P.D., Weaver, H.A., A'Hearn, M.F., Tozzi, G.-R, & Festou, M.C. 1996, 
DPS, 28, 09.29 

Meier et al. 1998, Science, 279, 842 

Millar, T.J., Bennet, A., & Herbst, E. 1989, ApJ, 340, 906 

Millar, T. J., Farquhar, P. R. A., & Willacy, K. 1997, A&AS, 121, 139 

Mumma, M.J. 1997, in "From Stardust to Planetesimals" , eds. Y. Pendleton & A.G.G.M. Tielens, 
(San Fransicso: Ast. Soc. of Pac. Conf. Series), 122, 369 

Neufeld, D.A. & Hollenbach, D.J. 1994, ApJ, 428, 170 

Pineau des Forets, G, Roueff, E., & Flower, D.R. 1989, MNRAS, 240, 167 

Reipurth, B., Bally, J., & Devine, D. 1997, AJ, 114, 2708 

Rodgers, S.D. & Millar, T.J. 1996, MNRAS, 280, 1046 



-9- 



Safronov, V.S. 1969, Evolution of the Protoplanetary Cloud and Formation of the Earth and 
Planets (Moscow: Nauka Press), [in Russian: english translation NASA TTF-667, 1972] 

Sandford, S.A. & Allamandola, L.J. 1990, Icarus, 87 

Schutte, W. 1998, in Laboratory Astrophysics and Space Research, eds. P. Ehrenfreund & H. 
Kochan (Dordrecht: Kluwer), in press 

Shalabiea, O.M., Casern, P., & Herbst, E. 1998, ApJ, 502, 652 

Sundstrom et al. 1994, Science, 263 , 785 

Talukdar, R.K. et al. 1996, J. Phys. Chem., 100, 3037 

Tielens, A.G.G.M. 1983, A&A, 119, 177 

Tielens, A.G.G.M. & Hagen, W. 1982, A&A, 114, 245 

Whittet, D.C.B. et al. 1998, ApJ, 498, L159 



This preprint was prepared with the AAS IAT^X macros v4.0. 



-10- 



Table 1. High Temperature Isotopic Reactions' 



Reaction 


A (cm 3 s x ) 


B 


C(K) 


Refs 


+ HD -> OD + H 


1.57 xl0~ 12 


1.7 


4639 


1 


+ HD -> OH + D 


9.01 xlO -13 


1.9 


3730 


1 


OH + HD -> H 2 + D 


2.12 xlO -13 


2.7 


1258 


2 


OH + HD -> HDO + H 


0.60 xHT 13 


1.9 


1258 


2 


OD + H 2 -» HDO + H 


1.55 xl0~ 12 


1.6 


1663 


3 



a Rate coefficient k = A(T/300 K) B e - c / T . 

1 Fit to theoretical predictions of Joseph, Truhlar, & Garret (1988). 
2 Talukdar et al. (1996) 

3 Rate same as OH + H 2 ; Talukdar et al. (1996). 
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Pre-Shock: T =30 K Shock: T =1000 K Post-Shock: T =30 K 
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Fig. 1. — Evolution of chemical abundances in our gas-grain three-stage model as a function of 
time. All abundances are relative to H2. This model includes molecular depletion and desorption 
from grain surfaces and the abundance of H2O, HDO, and CO2 on the grain surface is shown. Stage 
1 (first panel) represents the pre-shock stage with T gas = 30 K, Stage 2 (second panel) is the shock 
stage with a higher gas temperature (T gas = 1000 K; v s = 20 km s" 1 ), and Stage 3 (third panel) 
the post-shock stage with T gas = 30 K. The density and visual extinction is constant for all stages 
at n(H2) = 10 5 cm -3 and Ay = 10 mag. Note that the time axis for each panel covers a different 
range and shows the full evolution of chemical abundances from one stage to the next. 
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Fig. 2. — Plot of the grain surface [HDO]/[H20] against [C02]/[H20] shown as a function of time 
[log(t)] and gas and dust temperature (T). The calculated abundance ratios are shown only in the 
post-shock phase (stage 3). The points are the observed ratios in Comets Halley, Hyakutake, and 
Hale-Bopp. For Halley, Hyakutake, and Hale-Bopp the HDO data are from: Eberhardt et al. 
(1995), Balsiger et al. (1995), Bockelee-Morvan et al. (1998), and Meier et al. (1998), respectively. 
The CO2 measurements are from Kranowsky (1986), McPhate et al. (1996), and Crovisier et al. 
(1997). The box shows the range of representative values observed in the interstellar medium and 
[HDO]/[H20] ratios are from Gensheimer et al. (1996) and the [C02]/[H20] ratios in polar ices are 
from Schuttte (1998) and Whittet et al. (1998). 



